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Abstract 

In this paper, we study multivariate vector sampling expansions on general 
finitely generated shift-invariant subspaces. Necessary and sufficient condi- 
tions for a multivariate vector sampling theorem to hold are given. 
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1. Introduction and Main Results 

If H is a Hilbert space, we define =HxH x • • • x H(q term). Given 
/ = (/i,/2,-- • ,fq) T ,9 = (9i,92, ■ ■ ■ ,9q) T G , the inner product / and g 
is defined by 

q 

(f,9) H (i) = ^2(f P ,9p)H- 
P =l 

Let ifij = (<fij 7 i, fj,2, " ' > fj,r) T £ L 2 (R d )( r \ 1 < j < N be a stable generator 
for the shift-invariant subspace 

i.e., the sequence — a) : 1 < j < N, a € Z d } is a Riesz basis for V 2 . 

Recall that {<fij(- — a) ■ 1 < j < N,a G Z d } is a Riesz basis for V 2 , if there 
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exist two constants A, B > such that for any {c J)Q : 1 < j < N, a € Z rf } £ 
£ 2 (Z d )( N ) one has 
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3=1 a&L d 
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We assume throughout the paper that the vector functions in the shift- 
invariant subspace V 2 are continuous on W 1 . Equivalently(see [l[ or 0]), 
that the generator ipj , 1 < j < N is continuous on and 



iV 



sup x6Rd E E E \wA x - «)i 2 < °°- 

i=l p=i oez d 



If V^.(l < j < N) is defined by 



A? 



Then we have 

F 2 = V y 2 

i=i 

Define T^, : L 2 ([(j - 1)/N, j/N) d ) — > by 

= E Cp d^j{- 



a 



where c Fj , Q = iV d / 2 J^^/^ F(x)e 2 ^ T -'dz. 

For convenience, in this paper we make Xp( x ) = X[(p-i)/N,p/N) d ( x ) an d 
e Q (x) = N d / 2 e- 2nNiaT - x . 

Lemma 1.1. = EjLiT^ is an isomorphism between L 2 [0, l) d and V 2 . 

Proof. Since the sequence - a) : I < j < N,a eZ d ] is a Riesz 

basis for V 2 , Therefore, for any F 6 L 2 [0, l) d , we have 



It pii^ 
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B E \\ F ^lU-^/N,j/N] d 

i=i 
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L 2 [0,l) d - 
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Similarly, we also have 



\T„F 



> A\\F\ 



L 2 [0,l) d - 



□ 

Given a nonsingular matrix M with integer entries. Let 7^ + M Tr L d , 1 < 
k < m = |det(M)| be the m distinct elements of the coset space Z d /M T Z rf 
with 71 = 0. Define Q k = M- T ^ k /N + M~ T [0, l) d /N, 1 < k < m, we have 
(see 0, P. 110]) 



Q k n Q k , = and Vol ({J Q^j = i. 



Thus, for any function i 7 integrable in [0, 1/N) d and Z rf /iV-periodic, we have 

/[o,i/jv)«« = Er=i /q, ^oo^. 

Let ffj € L 2 [0, l) d , 1 < j < s, define 



G P (x) 



92{x)x P {x) 



and its related constants 
A G := 
Be : 



gi{x) Xp {x + M- T lm /N) 
g2(x) Xp {x + M- T lm /N) 

g s {x)x P (x + M- T lrn /N) 



,1 <p< N.(l.l) 



min essinf A min \G*Jx)GJx)] , 
max ess sup A max \GZ(x)G p (x)] . 

1<P<N x£[0jl/N )d 



Lemma 1.2. Suppose that gj € L 2 [0, l) d , 1 < j < s and G p (x), 1 < p < N 
is its associated matrix as in Then 

(a) The sequence |<fy- (x)x P (x)e a (M T x) : 1 < j < s,l < p < N, a £ Z rf j is 

a complete system for L 2 [0, l) d if and only if for any 1 < p < N the 
rank of the matrix G p {x) is m a.e. in [0, 1/N) d . 

(b) The sequence ^gj(x)x P (x)e a (M T x) : 1 < j < s,l < p < N, a £ Z d j is 

a bessel sequence for L 2 [0, l) d if and only if Bq < 00. In this case, the 
optimal Bessel bound is Bc/m. 

(c) The sequence ^gj(x)x P (x)e a (M T x) : 1 < j < s,l < p < N, a £ Z rf j is 

a frame for L 2 [0, l) d if and only if < Aq < Bq < 00. In this case, 
the optimal frame bounds is Ac/rn and Bc/m. 
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(d) The sequence ^gj(x)x P (x)e a (M T x) : 1 < j < s,l < p < N,a & Z d | is 
a Riesz basis for L 2 [0, l) d if and only if it is a frame and s = m. 

Proof. For any F £ L 2 [0, l) d , we have 
' F(x),g j (x)xp(x)e a (M T x 



L 2 [0,l) d 



[0,l) d 

F(x) gj (x)x P (x)e 

[o,i) d 

/ F(x)Xp(x)9j(x)x P (x) 
J[o,i) d 

[ F(x)xp(x)g J (x)xp(x)e^ maT - MTx dx 

J[0,l/N) d 



F(x)gj(x)x P (x) e a {M T x)dx 



2nNia T -M T x 



dx 



e 2nNia T -M T x^ x 



m „ 

= F ( x )Xp(x)gj(x)xp(x)e 
k=1 Q k 

= J 

J A 



2izNia T -M T x 



dx 



Y,(Fx P )(x + M~ T lk /N)x 

M-T[ ,l/N) d ^ 

(9jX P )(x + M- T lk /N)N d / 2 e^ NiaT - MTx dx. 
where we have considered the Z d /iV-periodic extension of F. Then 

s 

EE 



F(x),g j (x)x P (x)e a (M 1 



1 s 

= -Y, 



L 2 [0,l) d 
m 

^(F Xp )(x + M- T 7fc /iV)x 



k=i 



{gjXp){x + M~ T lk /N) 



L 2 (M- T [0,l/N) d ) 



(1.2) 



Denoting ¥ p (x) := ((Fx P )(x), • • • , (Fx P )(x + M~ T j m /N)) T , the above reads 



EE 



F(x),g j (x)xp(x)e a (M x) 



L 2 [0,l) d 



— \\G p (x)¥ p (x)\\ 2 . . M . 



(1.3) 
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Since 71, 72, • • • , 7 m are representatives of distinct cosets of Z d /M T Z d , there- 
fore for any 1 < k < m the matrix G(x + M~ T ^^/N) has the same columns 
of G(x). Hence, for any 1 < p < N rankG p (x) = m a.e. in [0, 1/N) d if and 
only if rankGp(x) = m a.e. in M~ T [0, l/N) d . Moreover, we have 

essinf X min \G*Jx)G p (x)] = essinf X min [G*(x)G p (x)] 

and 

ess sup X min [G* (x)G p (x)] = ess sup X min [G*(x)G p (x)] 

x€[0,l/N) d ' x£M- T [0,l/N) d 

To prove (a), assume that there exists set $7 C M~ T [0, 1/N) d with 
positive measure and 1 < po < N such that rankG Po (x) < m, x G O,. 
Then, there exists a measurable function v(x), such that G Po (x)v (x) = and 
\v(x)\ = 1 in a Define F £ L 2 [0, l) d such that ¥ Po (x) = v{x) if x G Q, 
¥ po (x) = if x € M~ T [0, l/iV) d \ft and F p (x) = if p ^ p . Hence, from 
(|1.3|) we obtain the system is not complete, conversely, if the system is not 
complete, by using (jl.3p we obtain a ¥ p (x) different from in a set with 
positive measure such that G p (x)¥ p (x) = 0. Thus rankGp(x) < m on a set 
with positive measure. 

If Bq < 00 then, for each F G L 2 [0, l) d , we have 

N s 

EEL 

p=l j=l 

JV 



F(x),5 i (j;)xp(x)e a (M T x) 



V — ||G p (x)F„(x)|| 2 , 

m 11 ; pv y|l L 2 (M- T [o,i/Ar) d ) 

p=i 



w 

A 

0) 

ivi - [u,i/iv )~ 1 

A? 



S G- E / E K^X* + M " T 7fc/^)| 2 

m p,JM-T [0 ,l/N) d ^ 
N rn 

S «-EE/ l(^X P )(x)| 2 dx 
m /n, 
p=l k=l ^ k 

N 

B G -E / \(F Xp )(x)\ 2 dx 
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,2 j_ Bq M ^,,2 



[(p-l)/N,p/N) d m 1 ' 



SC- 



Hence |g , j(x)xp(^)e Q ,(M T a;) : I < j < s,l < p < N,a £ Z d | is a bessel 
quence for L 2 [0, l) d and the optimal Bessel bound is less than or equal to 

m 

Let K < Bq. Then, there a set 9. K C M- T [0,l/N) d with positive 



measure and 1 < p < N such that A r 



G*Jx)G p (x) 



> K for x £ Q,k- 



Let F € L 2 [0, l) d such that its associated vector function ¥p(x) is if x € 

is an eigenvector of norm 1 associated with the 
largest eigenvalue of Gp(x)Gp(x) if x G VLk and ¥ p (x) = if p ^ p. Using 
(11,31). we obtain 



EE 



3=1 c 
> 



F(x), gj(x)xp(x)e a (M T x 



1 



L 2 [0,l) d 

K 



K\¥ § {x)\Ux = ^\\F\\l 2md . 



m JM- T [0,l/N) d 

Therefore if Bq = oo then ^gj{x)x p {x)e a (M T x) : 1 < j < s, 1 < p < 

N, a € Z d | is not a bessel sequence for L 2 [0, l) d , and if Bq < oo then the 

optimal Bessel bound is Bq/tti. This completes the proof of (b). The proofs 
of (c) are completely analogous. 

To prove (d), we assume that m = s and that the sequence is a frame. 
We see that it is a Riesz basis by proving that the analysis operator 

A : L 2 [0,l) d ^£ 2 {Z d )( Nxs \ 



i.e. 



A(F) :- 



F(x), gj(x)x P (x)e a (M x) 



L 2 [0,l) d 



1 <j < s,l <p < N,a e Z fl 



is surjective (see 0, Theorem 6.5.1]). To this end, notice that when m = s 
for any 1 < p < N the matrix G p , 1 < p < N is a square matrix and hence, 
the condition Aq > implies that for any 1 < p < N the inverse matrix 
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G p 1 (x) exists and its entries are essentially bounded. For 1 < p < N, let 
{c£a : 1 < j < s,a £ Z d } be an element of £ 2 (Z d )^ . For 1 < p < N, 1 < 
j < s we define the function 

a€Z d 

and let F be the function such that 

¥ p (x) = G p \x) (£[(*), ^(x), • • • ,e s (x)f ,x G M- T [0, l/iV) d . 

This function belongs to L 2 [0, l) d because the entries of G~ 1 (x) are essen- 
tially bounded. We have that G p (x)¥ p (x) = (^(x), £r? (x), ■ ■ ■ , £?(x)) T , and 
using (|1.2p we obtain that 



is a 



JM- T lO,l/N) d 

and consequently A(F) = |c^ a : 1 < j < s, 1 < p < iV, a G 

Conversely if {g j {x)x P {x)e a {M T x) : 1 < j < s, 1 < p < 2V, a G Z d | 

Riesz basis. Let {/j, p , a : 1 < j ' < s, 1 < p < 2V, a G Z d } be its dual Riesz 
basis. Then, by using (jl.2p we obtain for 1 < p < N 

,. m 

mSafiSjji = / ^2(f j ,, pfi x P )(x + M~ T ~f k /N) x 

■/M- T [0,l/AO d ^ 

Therefore, for 1 < j, j' < s, we have 
m 

^(/j', P ,oX P )(x + M- T 7fe /^)( 5j Xp)(x + M~ T lk /N) = md j>f ,a.e. 

k=l 

Thus the matrix G(x) has a right inverse; in particular, s < m. As a 
consequence (a) we have s > m and, finally, s = m. □ 
We consider s linear- invariant systems Cj,l < j < s in L 2 (M d )^ such 
that for any / = (/i, / 2 , • • • , / r ) T G L 2 (M d )M, 



(£,•/)(*) = [/ * P](i) = W f q {x)p hq {t - x)dx, 

U jRd 
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where P(x) is an sxr matrix with entries pj :Q S L 1 (M. d ), 1 < j < s, 1 < q < r. 

T 



Let hj(t) = \p jtl {-t),pj^{-t), ■ ■ ■ ,pj t r(-t)\ , we have 
(C 3 f)(t) = {f(-),h J (.-t)) LHRd)(r) . 

The set of systems {C\, £2, • • • , £ s } is an M-stable filtering sampler for 
if there exist two positive constants C\ and C2 such that [H] for any 
/ = /(i) + /( 2 ) + . . . + fW G where G Vg,, we have 



p=i j=i aez d 



L 2 (K d )('')- 



For 1 < j < s and 1 < p < 2V, we define gj )P {x) by 

: = E (^ P )(»)e- 2 - NiaTx . (1.4) 

Lemma 1.3. Xei / 6e a function in V£ such that f = /W + /^H h/^ 

w/iere € V* p and = T Vp F p ,F p € L 2 ([(p - 1)/N,p/N) d ). For every 
1 < j < S; we have 

(Cjf^)(Mp) = (F p (.),^)e^(M T -)) L2[oi)d (1.5) 

Proof. For each £ Z d we have 
(^•/W)(M/3) = (/W(.),/ ij (--M/3)} L2(Rd)M 

= ( E c Fp^v(- ~ a )i h j (• - M P) ) 

\a& d I T,2(md\(r) 



E °F P , a {<P P {- - a),hj (• - Mj3)) L2(Rd)(r) 

a€Z d 

E o Fp , a {C jVp ){Mp-a) 

E ^p(')' ea ('))L 2 ([(p-l)/Ar,p/AT) d 

} (/: iV 9 p )(M/3-a) 

aGZ d 

( E T^Pp){M0 - a)e a (-)\ 

\ ' L 2 ([(p-1)/Af,p/Af) d ) 

□ 
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Theorem 1.4. Assume that the function gj tP (x) given in ji.^p belong to 
L°° ([(p - l)/N,p/N) d ) /or eac/i 1 < j < s and 1 < p < N. Let G p (x) be 
the associated matrix define in [(p — 1)/N,p/N) d as in The following 

statements are equivalents: 

(a) Aq > 0; 

(b) The set of systems {Ci, £2, ■ •■ , As} is an M -stable filtering sampler 
for V%; 

(c) For 1 < p < N, there exist vectors (d^(x), d\{x), ■ ■ ■ , d%{x)) with en- 
tries d P - G L°° ([(p - 1)/N,p/N) d ) satisfying 

«(x), dP 2 (x), ■ ■ • , d p (x))G p (x) = (1, 0, • • • , 0) 

a.e.in[(p-l)/N,p/N) d ; (1.6) 

(d) There exists a frame for having the form {Sj(t — Ma) : 1 < j < 
s, 1 < p < N, a G Z d } such that for any f G V<% 

s N 

f = m ^2Yl Yl ^j&\Ma)^{t -Ma) mL 2 (M d )W. (1.7) 
j=\ p =i ae %d 

Proof. Part (c) in Lemma 11.21 proves that conditions (a) and (b) are 
equivalent. 

If Aq > then for any 1 < p < N, ess inf X £\o t i/N) d det [G*(x)G p (x)] > 
and consequently, there exists the pseudo-inverse matrix 

gI(x) = [g;(x)g p ( x )]- 1 g;(x). 

Moreover, its entries are essentially bounded and its first row satisfies (jl.6p . 
Therefore, (a) implies (c). 

Next, we will prove that the condition (c) implies (d). Since we have 
assumed that gj, P (x) G L°° ([(p — 1)/N,p/N) d ) for any 1 < j < s and 
1 < p < N, Lemma ll.2( b) proves that 

{g jtP (x)e a (M T x) : 1 < j < s, 1 < p < N, a G Z d ] 

is a Bessel sequence in L 2 [0, l) d . The same argument proves that 

{md P (x)e a (M T x) : 1 < j < s, 1 < p < N, a G Z d } 
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is also a Bessel sequence in L 2 [0,l) d . By CCU) and CCS]), these two Bessel 
sequences satisfy 



TV 



F(x) = mY,J2 E (F(-),9jA-)^(M T -))dF j (x)e a (M T x),F e L 2 [0,l) d . 

3=1 p=l aGZ d 

Hence, they form a pair of dual frames for L 2 [0, l) rf (see[3, Lemma 5.6.2]). 
Since Sj(t — Ma) = T^[<i^(-)e a (M T -)](t) and is an isomorphism, the 
sequence - Ma) : l<j<s,l<P<N,aeZ d } is a frame for V^. 

Last, we prove that the condition (d) implies (6). Notice that since 

we have assumed that ^gj tP (x)e a (M T x) : 1 < j < s,l < p < N, a € Z rf j 
a Bessel sequence with bound Bq/tti and 



is 



{L ] f^){Mp) = {F{.),g hp {-)e p {M 1 '-) 



L*([(p-1)/N,p/N)*) 



For each / e VJ?, we have 



TV s 
P=l i=l agZ d 



^ — l^llL2[0,l) d ^ 



^gIIt- 1 



1 1 oper n n 2 



III 



If {5j(t - Ma) : 1 < j < s, 1 < p < N, a € Z d } is a frame for V^, then the 
formula (|1.7p gives 



L 2 (R d )M 



m 



s TV 

EEE £ ^ (P) ( M «) 5 ^- M «) 



L 2 (R d )M 



TV s 

< - 2 ^EEE|^/ (P) ( M ^ 

p=i i=i oez d 

where C is a Bessel bound for {S?(t - Ma) : 1 < j < s, 1 <p < N,a € Z d }. 
Hence, the set {C\, C2, • • • , £ s } is an M-stable filtering sampler for V 2 . □ 
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